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Using  in  situ  high-resolution  synchrotron  x-ray  scattering,  the  Curie  temperature  has  been  deter¬ 
mined  for  ultrathin  c-axis  epitaxial  PbTi03  hlms  on  conducting  substrates  (SrRu03  on  SrTi03),  with 
surfaces  exposed  to  a  controlled  vapor  environment.  The  suppression  of  Tc  was  relatively  small,  even  for 
the  thinnest  film  (1.2  nm).  We  observe  that  180°  stripe  domains  do  not  form,  indicating  that  the  de¬ 
polarizing  field  is  compensated  by  free  charge  at  both  interfaces.  This  is  conhrmed  by  ab  initio  calcu¬ 
lations  that  find  polar  ground  states  in  the  presence  of  ionic  adsorbates. 
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One  of  the  outstanding  fundamental  questions  for  ferro¬ 
electric  thin  films  is  understanding  the  stability  of  the  polar 
phase  when  the  polarization  has  a  component  perpendicu¬ 
lar  to  the  film  plane.  Although  such  orientations  are  desired 
for  many  applications  [1]  and  can  he  obtained  (e.g.,  by 
compressive  epitaxial  strain  [2]),  the  polarization  change  at 
the  film  boundaries  creates  a  “depolarizing  field”  that 
must  be  neutralized  for  the  polar  phase  to  be  stable  [3]. 
Two  mechanisms  are  available  to  reduce  the  depolarizing 
field  energy:  compensation  by  free  charge  at  the  bounda¬ 
ries  [4]  or  the  formation  of  equilibrium  stripe  domains  with 
oppositely  oriented  polarization  [5].  In  both  cases,  the 
trade-off  between  bulk  energy  gain  and  surface  energy 
cost  leads  to  a  suppression  of  the  phase  transition  to  the 
polar  phase  as  films  become  thinner.  These  fundamental 
size  effects  may  dramatically  alter  behavior  in  ultrathin 
films. 

Here  we  focus  on  stabilization  of  the  single-domain  state 
in  ultrathin  ferroelectric  perovskite  films  by  interfacial 
charge.  Key  issues  are  determining  the  critical  thickness, 
below  which  the  monodomain,  perpendicularly  polarized 
state  is  not  stable,  and  understanding  the  nature  of  the 
interfacial  charge  compensation.  Previous  experimental 
studies  of  ultrathin  epitaxial  films  on  conductive  substrates 
have  been  carried  out  using  electrical  measurements  [6], 
piezoresponse  force  microscopy  (PEM)  [7,8],  and  x-ray 
lattice  parameter  measurements  [8].  These  studies  have 
found  stable  monodomain  ferroelectricity  at  room  tem¬ 
perature  in  films  as  thin  as  3-4  nm.  Lattice  parameter 
measurements  indicated  a  50%  reduction  in  PbTi03  polar¬ 
ization  when  film  thickness  was  reduced  to  3  nm  [8],  which 
was  explained  using  a  model  based  on  the  finite  screening 
length  for  charge  in  conducting  electrodes  [4].  Interest¬ 
ingly,  in  some  of  these  studies  [7,8],  the  film  did  not  have 
an  electrode  on  the  exposed  surface.  This  suggests  that 
interface  compensation  adequate  to  stabilize  ultrathin  polar 


films  can  occur  by  a  mechanism  other  than  electronic 
conduction,  e.g.,  by  the  accumulation  of  charged  ions 
[9].  Eor  relatively  thick  ferroelectric  films  exposed  to 
ambient  atmosphere,  there  is  strong  experimental  evidence 
for  surface  compensation  by  ionic  adsorption  [7,10-12]. 
Eurthermore,  monodomain  ferroelectric  films  have  been 
observed  on  nonconducting  substrates  [13-15],  and  it  has 
been  proposed  that  high-mobility  electrons  can  exist  at 
heterointerfaces  in  perovskite  insulators  due  to  unusual 
bonding  states  [16].  The  adequacy  of  these  alternative, 
“chemical”  rather  than  “electronic”  mechanisms  for  re¬ 
ducing  depolarization  field  energy  sufficiently  to  stabilize 
ultrathin  polar  films  is  relatively  unexplored. 

In  this  Letter,  we  report  a  joint  experimental  and  theo¬ 
retical  study  of  the  paraelectric -ferroelectric  phase  transi¬ 
tion  in  ultrathin  PbTi03  films  grown  on  conducting 
substrates,  with  surfaces  exposed  to  a  controlled  vapor 
environment.  The  structures  of  films  as  thin  as  1.2  nm 
(3  unit  cells)  were  determined  as  a  function  of  temperature 
(T)  using  high-resolution  grazing-incidence  x-ray  scatter¬ 
ing.  The  polar  phase  was  observed  to  form  and  remain 
single  domain  for  all  thicknesses  studied  (1. 2-9.2  nm). 
The  observed  stability  of  the  monodomain  phase  is  ex¬ 
plained  by  ab  initio  calculations  that  find  a  polar  ground 
state  if  ionic  adsorbates  are  present  on  the  surface.  The 
polarization  direction  is  predicted  to  depend  on  the  chemi¬ 
cal  nature  of  the  adsorbate. 

The  conducting  substrates  were  epitaxial  SrRu03  films 
grown  on  SrTiO3(001)  single  crystals  by  pulsed  laser 
deposition  [17].  The  PbTi03  films  were  grown  by  metal- 
organic  chemical  vapor  deposition  (MOCVD)  under  con¬ 
ditions  described  previously  [18].  Grazing-incidence  syn¬ 
chrotron  x-ray  scattering  was  used  for  in  situ  monitoring  of 
the  film  during  epitaxial  growth  and  for  subsequent  obser¬ 
vation  of  the  phase  transition  as  a  function  of  T  [14].  This 
method  allows  study  at  high  T  while  maintaining  film 
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Stoichiometry  and  a  well -controlled  vapor  interface  and  the 
study  of  ferroelectricity  in  films  too  thin  to  he  characterized 
hy  other  methods.  Eour  PbTi03  films  of  different  thickness 
were  studied:  Three  of  the  films  (1.2, 2.0,  and  3.6  nm)  were 
grown  at  930  K  on  10-nm-thick  SrRu03  layers;  the  thickest 
PbTi03  film  (9.2  nm)  was  grown  at  990  K  on  a  50-nm-thick 
SrRu03  layer.  All  of  the  layers  replicated  the  high  crystal¬ 
line  quality  of  the  substrates  (0.01°  typical  mosaic)  and 
were  measured  to  be  fully  lattice  matched  to  the  underlying 
SrTi03.  This  compressive  epitaxial  strain  produces  polar¬ 
ization  perpendicular  to  the  film  plane  [2]. 

X-ray  experiments  were  carried  out  at  beam  line  12ID-D 
of  the  Advanced  Photon  Source,  as  described  previously 
[14].  The  Curie  temperature  was  determined  by  mea¬ 
suring  the  T  dependence  of  the  c  (out-of-plane)  lattice 
parameter.  Eigure  1  shows  measurements  along  the  30L 
crystal  truncation  rod  (CTR)  for  a  9.2  nm-  (23  unit-cell-) 
thick  PbTi03  film  at  various  temperatures.  This  region  of 
the  CTR  extends  through  the  304  peaks  of  PbTi03, 
SrRu03,  and  SrTi03.  Erom  the  position  of  the  PbTi03 
peak,  one  can  see  that  c  increases  as  the  film  polarizes 
below  7(3,  because  of  the  strong  polarization-strain  cou¬ 
pling  in  this  system  [2].  The  many  fringes  observed  from 
the  thicknesses  of  the  SrRu03  and  PbTi03  layers  indicate 
the  high  quality  of  the  interfaces. 

Eor  the  9.2  nm  film,  we  also  performed  ex  situ  room-7 
PEM.  We  were  able  to  “write”  inverted  domains  by  apply¬ 
ing  a  positive  voltage  to  the  tip  and  to  determine  that  the  as- 
grown  state  of  the  film  was  a  single  domain  having  a 
polarization  directed  out  of  the  film  (“up”  polarization). 
Eor  all  films,  we  searched  in  situ  near  for  satellites  in  the 
diffuse  x-ray  scattering  around  the  PbTi03  Bragg  peaks, 
which  occur  when  equilibrium  180°  stripe  domains  are 
present  [5,14].  None  were  observed,  indicating  that  all  the 
films  transform  directly  into  the  monodomain  polar  state. 


FIG.  1.  X-ray  scattering  intensity  along  the  307  CTR,  showing 
the  304  Bragg  peaks  and  finite-thickness  fringes  from  9.2-nm- 
thick  PbTi03  and  50-nm-thick  SrRu03  films  at  various  tempera¬ 
tures.  The  shift  of  the  PbTi03  peak  to  lower  upon  cooling 
indicates  an  expansion  of  the  c  lattice  parameter. 


Values  for  the  PbTi03  lattice  parameter  c  as  a  function 
of  7  were  determined  by  fitting  the  x-ray  CTR  data  to  a 
3-layer  model  (two  films  and  a  substrate)  [13].  Eor  the  fits 
shown  here,  we  assumed  that  the  PbTi03  polarization 
direction  is  up,  with  the  square  of  atom  displacements 
proportional  to  the  change  in  c  from  its  7  =  0  value,  and 
that  the  SrRu03  is  SrO  terminated  [19]  and  the  PbTi03  is 
PbO  terminated  [20].  Changing  these  assumptions  does  not 
significantly  affect  the  results  for  c.  Eor  each  sample,  the 
numbers  of  PbTi03  and  SrRu03  unit  cells  were  fixed. 
Eight  parameters  were  varied  in  the  fitting  procedure: 
lattice  constants,  layer  roughnesses,  and  interface  offsets 
for  the  SrRu03  and  PbTi03  layers,  a  scale  factor  for  the 
SrRu03,  and  an  overall  scale  factor  for  the  total  scattered 
intensity.  This  simple  structural  model  is  able  to  reproduce 
the  CTR  intensities  very  well  for  all  thicknesses  and  tem¬ 
peratures.  Typical  best  fits  are  shown  with  the  data  in 
Pig.  1. 

Eigure  2  shows  the  dependence  of  c  on  7  for  each  film 
thickness,  as  well  as  the  predicted  7  dependence  of  c  for 
thick,  coherently  strained  PbTi03  lattice -matched  to 
SrTi03  [2,14].  In  the  thick  limit,  the  (second  order)  non¬ 
polar  to  polar  transition  is  located  at  7“  =  1025  K,  as 
indicated  by  the  abrupt  change  in  slope.  We  extracted  7^ 
for  each  film  by  estimating  the  temperature  at  which  the 
slope  of  c(7)  changed,  as  shown  by  the  solid  lines  in  Pig.  2. 
The  theoretical  slope  was  used  above  T(^.  Since  we  have  no 
data  above  7(3  for  the  3.6  nm  film,  that  estimate  of  7(3  is  a 
lower  bound. 

These  results  stand  in  some  contrast  to  recently  reported 
room-7  lattice  parameters  of  PbTi03  films  on  Nb-doped 
SrTi03  substrates  [8],  which  showed  a  decrease  in  c  for 
thinner  films.  We  see  no  change  in  c  as  a  function  of 
thickness  at  our  lowest  7  (550  K),  even  for  smaller  thick¬ 
nesses  than  previously  measured,  although  there  is  a  clear 
variation  in  7(3  with  thickness.  Note  also  that  we  observe  a 
systematic  change  in  c  with  thickness  in  the  nonpolar 


FIG.  2.  The  c  lattice  parameter  vs  7  for  each  film,  with  solid 
lines  showing  the  break  in  slope  used  to  estimate  Tq.  The  upper 
curve  is  the  Landau  theory  prediction  for  thick  PbTi03  on 
SrTi03  [2]. 
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phase.  Caution  should  thus  he  used  when  relating  c  to 
polarization  or  suppression  in  ultrathin  films  using 
measurements  at  only  one  T. 

The  values  of  as  a  function  of  film  thickness  are 
shown  in  Eig.  3.  Also  shown  are  values  of  for  PbTi03 
grown  directly  on  insulating  SrTi03,  in  which  180°  stripe 
domains  have  formed  [5].  In  both  cases,  increases 

towards  the  calculated  7“  for  thick  films  and  decreases 
by  hundreds  of  degrees  below  this  limit  for  ultrathin  films. 
The  transition  temperatures  for  PbTi03  on  SrRu03  are 
somewhat  higher  than  those  on  SrTi03,  in  agreement 
with  the  observed  difference  in  the  equilibrium  domain 
morphology:  The  conducting  SrRu03  electrode  lowers  the 
energy  of  the  single-domain  state,  producing  a  direct  tran¬ 
sition  from  unpolarized  to  monodomain  polarized  at  a 
higher  7^.  The  dependence  of  7“  —  7^  on  film  thickness 
we  observe  does  not  obey  the  power  law  predicted  by  the 
simplest  theories  of  screening  in  electrodes  [21]  but,  in¬ 
stead,  appears  to  be  similar  to  that  previously  seen  for  films 
on  insulating  substrates. 

The  observation  that  polarization  occurs  without  forma¬ 
tion  of  180°  domains  implies  that  both  the  top  and  bottom 
interfaces  of  these  films  must  be  almost  completely  com¬ 
pensated  by  free  charge  [3].  Thus,  a  SrRu03  electrode  will 
compensate  PbTi03  with  a  voltage  offset  small  enough  that 
even  a  3-unit-cell-thick  film  polarizes  above  room  tem¬ 
perature.  This  is  in  agreement  with  recent  ab  initio  calcu¬ 
lations  on  films  with  two  SrRu03  electrodes  [22].  The 
behavior  of  the  top  surface  is  perhaps  even  more  intrigu¬ 
ing — although  there  is  no  conductor  to  supply  electronic 
charge,  exposure  to  the  ambient  vapor  of  the  MOCVD 
growth  environment  evidently  supplies  sufficient  free 
charge  from  ions  to  neutralize  the  depolarizing  field. 

To  quantitatively  verify  whether  adsorbed  ions  can  sta¬ 
bilize  the  polar  state  in  ultrathin  films,  we  have  performed 
ab  initio  density  functional  theory  (DET)  calculations  on 
the  structure  and  energetics  of  PbTi03  films  on  SrRu03 
bottom  electrodes  with  various  molecules  adsorbed  to  the 
surface.  Based  on  the  composition  of  the  MOCVD  envi- 
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ronment,  we  chose  to  study  H,  O,  OH,  H2O,  and  CO2 
adsorption.  Calculations  were  performed  using  methods 
described  previously  [22,23].  We  modeled  3-unit-cell- 
(L2-nm-)  thick  PbTi03  films  supported  by  3  unit  cells  of 
SrRu03,  with  a  SrO-Ti02  interface  between  PbTi03  and 
SrRu03  and  a  PbO  terminated  PbTi03  surface.  A  vacuum 
of  >2  nm  separated  periodic  copies  of  the  structures  in  the 
direction  perpendicular  to  the  surface,  and  a  dipole  correc¬ 
tion  was  included  in  the  center  of  the  vacuum  region  to 
remove  the  artificial  electric  field  due  to  the  asymmetry  in 
“electrode”  materials.  The  in-plane  lattice  constant  was 
fixed  to  the  calculated  zero-stress  a  value  for  bulk  PbTi03, 
which  approximates  the  effect  of  the  SrTi03  substrate  in 
the  experiments.  Data  are  reported  for  films  with  one 
adsorbate  on  a  single  unit-cell  surface  area;  modeling  of 
the  c(2  X  2)  surface  showed  that  the  experimentally  and 
theoretically  determined  reconstruction  [20,24]  was 
present  both  with  and  without  adsorbates  but  did  not  sig¬ 
nificantly  affect  the  energy  difference  between  the  bare  and 
adsorbed  states. 

Eor  a  bare  surface,  the  ground  state  of  the  ferroelectric 
film  was  found  to  be  nonpolar.  With  OH,  O,  or  H  adsor¬ 
bates  present,  however,  the  ground  state  was  polar,  with 
atoms  displaced  from  the  centrosymmetric  positions. 
Table  I  gives  values  of  polarization  relative  to  the  theoreti¬ 
cal  bulk  value  Pbuik  =  0.75  C/m^,  calculated  from  the 
average  over  all  PbO  and  Ti02  layers  of  the  cation-anion 
displacements  Az  =  Zcation  “  ^o>  s^^h  divided  by  the  cor¬ 
responding  Az  in  bulk  theoretical  PbTi03.  Since  the  con¬ 
ductive  SrRu03  electrode  can  provide  either  positive  or 
negative  compensation  charge,  the  chemical  nature  of  the 
adsorbate  determines  the  direction  of  polarization.  An 
overlayer  of  OH  or  O,  which  bind  to  the  surface  Pb, 
enforces  an  upwards  polarization,  while  an  overlayer  of 
H,  which  bind  to  the  surface  O,  stabilizes  polarization  in 
the  opposite  direction.  On  the  other  hand,  CO2  adsorption 
gives  a  very  weak  polarization,  and  undissociated  H2O 
molecules  bind  only  weakly  to  the  surface,  preserving 
the  nonpolar  state. 

Table  I  also  shows  the  differences  in  energy  AEopi- 
between  the  film  with  bound  adsorbate  and  the  separated 
bare  film  and  free  adsorbate  atom  (or  OH  molecule).  To 
compare  these  7  =  0  K  energies  to  experimental  condi¬ 
tions,  we  estimated  the  standard  Gibbs  free  energy  of 
adsorption  AG°(7)  =  AH^(T)  —  TAS^iT).  The  change 
in  enthalpy  AH^  is  estimated  from  AE’^pj  by  adding 
PAV  ~  —kgT.  A  small  correction  for  zero-point  energy 


TABLE  I.  Polarization  and  AZt^pr  from  DFT  calculations, 
with  estimated  reaction  energies  per  adsorbate  at  7  =  300  K. 


Adsorbate 

7/|7buikl 

AZ?dft 

(eV) 

AH° 

(eV) 

AG° 

(eV) 

AG 

(eV) 

OH 

0.7 

-2.00 

-2.34 

-1.96 

-0.20 

0 

0.9 

-0.92 

-0.95 

-0.52 

1.95 

H 

-0.8 

-2.90 

-2.93 

-2.65 

0.79 
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and  spin  polarization  is  also  applied  for  the  case  of  OH 
adsorption.  The  change  in  entropy  A  5*^  is  primarily  due  to 
the  differences  between  the  hound  and  free  adsorbates.  We 
estimated  the  entropy  of  the  bound  adsorbates  [25]  and 
used  tabulated  entropies  for  free  OH,  O,  and  H  at  1  bar 
[26].  Values  of  AG®  at  T  =  300  K  are  given  in  Table  I.  The 
free  energy  change  under  experimental  conditions  is 
AG  =  AG®  —  where  Pexp  is  th®  adsorbate  par¬ 

tial  pressure.  Eor  these  experiments,  the  partial  pressure  of 
O2  was  controlled  at  3.3  X  10“^  bar,  while  that  of  H2O 
varied  between  2  X  10“^  and  1  X  10“®  bar  depending 
upon  the  reactions  of  the  MOCVD  process  gases.  To  obtain 
the  values  of  AG  in  Table  I,  the  partial  pressures  of  free  O, 
H,  and  OH  at  T  =  300  K  were  calculated  using  PhjO  ^ 
5  X  10“^  bar  [26].  Inspection  of  these  AG  values  shows 
that  only  OH  adsorption  is  thermodynamically  favored  in 
our  experimental  environment. 

The  thermodynamic  stability  of  OH  adsorbates  implies 
that  the  stable  monodomain  polarization  direction  is  up,  in 
agreement  with  the  PEM  result.  The  simple  thermody¬ 
namic  estimate  presented  here  suggests  that  the  full  cover¬ 
age  of  OH  would  desorb  from  the  3-unit-cell-thick  film 
above  T  =  450  K,  resulting  in  a  nonpolar  film.  Prelimi¬ 
nary  calculations  that  allow  for  partial  coverage  of  OH  give 
a  higher  transition  T,  in  closer  agreement  with  the  observed 
Tc  =  700  K. 

Our  DPT  calculations  verify  that  surface  charge  passi¬ 
vation  by  OH  adsorbates  is  indeed  adequate  to  stabilize  the 
observed  monodomain  state  in  these  films.  Eor  ultrathin 
films,  the  chemisorption  energy  of  2  eV  is  much  larger  than 
the  bulk  free  energy  difference  between  the  polar  and 
nonpolar  states  at  zero  field  (e.g.,  0.2  eV  per  unif-cell 
area  for  a  3-unil-cell-lhick  film).  Such  a  sfrong  influence 
of  inlerfacial  chemisfry  on  polarizafion  has  broad  implica- 
fions  for  ulfrafhin  ferroeleclric  films.  To  undersfand  polar- 
izafion  sfabilify,  if  will  be  necessary  for  measuremenfs  fo 
be  performed  wifh  confrolled  inlerfacial  chemisfry.  The 
behavior  of  ferroeleclric  films  wifh  exposed  surfaces 
[7,8,10-12]  may  differ  significanlly  from  Ihose  sand¬ 
wiched  belween  Iwo  eleclrodes  [4,6,22].  In  Ihe  former 
case  (e.g.,  PPM  experimenls),  Ihe  swilching  mechanism 
likely  involves  a  change  in  adsorbate.  Charged  impurilies 
al  buried  interfaces  may  play  a  similar  chemical  role  lo  Ihe 
surface-adsorbed  ions  considered  here.  Exploilalion  of 
Ihese  effecls  for  novel  devices  or  lemplaling  techniques 
may  be  possible,  bolh  Ihrough  chemical  conlrol  of  polarity 
and  Ihrough  polarizafion  conlrol  of  ionic  adsorplion. 
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